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<^ ■ ABSTRACT 
< 



We show that X-ray clusters would have cooled substantially over a Hubble 
time by transport of heat from their hot interior to the their envelope, if the heat 
conductivity had not been heavily suppressed relative to the Spitzer value due to 



thermal conduction can not be the mechanism that prevents cooling flows there. 



^ ■ magnetic fields. The suppression is required in order for the observed abundance 

^ ■ of hot X-ray clusters to be consistent with predictions from popular cosmological 

Tij- . models. If a similar or stronger suppression factor applies to cluster cores, then 
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■ 1. Introduction 

Recently, the old idea that heat conduction may suppress cooling flows in X-ray clusters 
^ \ (Binney & Cowie 1981; Tucker & Rosner 1983; Bertschinger & Meiksin 1986; Bregman & 

.c^ ; David 1988; Gaetz 1989; Rosner & Tucker 1989; David et al. 1992; Pistinner & Shaviv 
1996; Dos Santos 2001) has been revived due to the apparent lack of strong cooling flows 
in Chandra and XMM-Newton data (Fabian et al. 2001; Bohringer et al. 2001; Molendi & 
Pizolato 2001). It was argued that if the heat conductivity is not suppressed by more than 
a factor of a few relative to the Spitzer value, then the inward heat flow due to the positive 
temperature gradient in cluster cores would be sufficient to compensate for the energy loss 
caused by Bremsstrahlung cooling of the gas (Narayan & Medvedev 2001; Gruzinov 2002; 
Voigt et al. 2002). 

The purpose of this paper is to caution that a large heat conduction coefficient would 
also lead to dramatic cooling of the entire cluster gas due to energy transport outwards into 
the cooler, surrounding intergalactic medium. By analyzing ASCA data for 30 X-ray clusters, 
Markevitch et al. (1998) have identified a composite temperature profile that declines by a 
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factor of ~ 2 from the core out to about half the virial radius (see also Finoguenov et al. 
2001). This profile has been confirmed in a recent BeppoSAX data set of 21 clusters (De 
Grandi & Molendi 2002). Here we use this profile to calculate the resulting conductive heat 
transfer from the cluster interior to the surrounding envelope. 

We note that the above ASCA results were challenged in several clusters by preliminary 
XMM-Newton data which showed a nearly isothermal profile (Arnaud et al. 2001a,b; Pratt 
et al. 2001). However, since clusters represent the hottest structures in the universe there 
is no doubt that they are surrounded by cooler gas at a sufficiently large radius. In fact, 
if the temperature profile is nearly fiat all the way out to the radius where the surface 
brightness of the cluster is too faint to be detectable, then the temperature gradient around 
the (more distant) virial radius of the cluster must be even steeper than assumed here (since 
the temperature must eventually approach the ambient value over a shorter range of radii) . 
Hence, the heat flux that we calculate based on the Markevitch et al. (1998) temperature 
proflle provides a conservative lower limit for the conductive energy loss of cluster interiors. 

Infall of gas onto the cluster could in principle suppress conduction of heat outwards. 
However, numerical simulation show that much of this infall is conflned within fllaments 
that cover only a small fraction of the surface area of the accreting object, and the gas flow 

is often clumpy and episodic (e.g., Yoshida et al. 2002, and references therein). Hence, 
conduction is unlikely to be suppressed by infall in all directions at all times. In this paper 
we focus our attention on the virialized region of the cluster, where the gas is stationary 
and in hydrostatic equilibrium. Since most of the mass is located in the outer envelope of 
the cluster, a considerable transport of heat out of the cluster core would only have a minor 
effect on the cluster envelope. For example, a reduction by ~ 50% in the X-ray temperature 
of the cluster - which is dictated by the core at ~ 300 kpc - would typically only lead to 
a ~ 5-10% increase in the outer envelope temperature at a radius of ~ 3 Mpc. In the 
following we examine the question of whether the cores of clusters, which dominate their 
X-ray emission, would cool through thermal conduction by more than they are allowed to 
cool over a Hubble time, given the observed temperature profiles of Markevitch et al. (1998). 
Throughout the discussion we assume a fully-ionized, hydrogen-helium plasma with a helium 
mass fraction of 24%, for which the mean particle mass (including the electrons) is /j, = 0.59 
in units of the proton mass rrip. 

2. Cooling of X-ray Clusters by Heat Conduction 

The time over which thermal energy is drained out of an X-ray cluster by conduction 
can be found by dividing the total thermal energy of the cluster, (Mg//imp)(|A;sr), by the 
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rate of conductive heat loss across its boundary radius r^, 



cond — 



\{Mg/^imp)kBf 
^7:Tl[K\dT/dr\V 



where T = M^^ JJ*" T{r)pg{r)4:7ir^dr is the (mass- weighted) mean temperature of the cluster, 
Pg{r) is the mass density of the gas, Mg = JJ*" Pg{r) Anr^dr is the total gas mass out to the 
radius r^, k. is the coefficient of heat conductivity, and ks is Boltzmann's constant. The 
hydrostatic equilibrium equation, (GM^o^pg/r^) — —dripghsT / iinip) , yields 



^ GuTTLpJ \^ 91nr 91nr 



dlnpg ^ dlnT 



(2) 



where fg = (Mg/Mtot) is the gas mass fraction at r^, G is Newton's constant, and dr = 
By substituting equation (2) into equation (1) we get 
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Thus, the characteristic coohng time for the entire cluster depends on the boundary values 
of the conductivity coefficient, ^(r;,), and the effective "adiabatic index" of the gas, 7efr = 

[l + {dr\nT/dr\npg)]r,. 

As customary in the literature, we identify the boundary of the hot interior of the 
cluster as the radius, rigo, where the average interior density of the gas is 180 times the 
mean cosmological density. The universal temperature profile derived by Markevitch et 
al. (1998) from ASCA data implies [91nT/(91nr]^, ^ -0.4, = [T{n)/f] ^ 0.6, and 
7cfr ~ 1.24. In order to get a numerical value for the cooling time tcond, we normalize the 
conductivity coefficient by the Spitzer (1962) value, nsp{^b) = 5 x 10^^/cB(/3f,Tio)^/^ cm~^ s~^ 
where Tio = (fcsT/lO keV), and define the suppression fractor r] = k/kq^. By substituting 
the above parameter values into equation (3) we get^, 

-3/2 



cond 



knT 



Gyr. 



10 keV, 

Over a cluster lifetime, Td, the cluster temperature is expected to decline by a fraction 



AT 



^cond 



0.2 



0.1 



10 Gyr) VO.15 



kef 
10 keV 



3/2 



(4) 



(5) 



-'^Note that the heat flux saturates when the Coulomb mean- free-path Acoui becomes longer than the 
characteristic scale of the temperature variation [dlnT / dr)~^ (Sarazin 1988). For a conductivity coefficient 
which is lower than the Spitzer value by a factor of r], the effective mean- free-path (at the same thermal speed) 
is AeflF ~ r?Acoui- From the T-vs-riso relation of Evrard et al. (1996) we get (Aeflp/rigo) ~ 0-l(^/0-l)^fo^- 
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where we have assumed that Td <S ^cond- Note that the cluster hfetime cannot be shorter 
than the sound crossing time across its diameter, ~ 6(ri8o/5 Mpc)!"^^^^ Gyr. Coohng is 
inevitable (and cannot be eliminated by adiabatic heating of the gas as it settles towards 
the center) since the gravitational potential that confines the gas is dominated by the dark 
matter. 

Equation (5) implies that hot X-ray clusters in the local universe must have been even 
hotter when they formed. The abundance of X-ray clusters is exponentially suppressed at 
high temperatures (Henry & Arnaud 1991; Eke ct al. 1996; Viana & Liddle 1996; Pierpaoli 
et al. 2001; Ikcbc ct al. 2002) and even more so at earlier cosmic times (Fan, Bahcall, & Cen 
1997; Bahcall & Fan 1998; Evrard et al. 2002). It therefore becomes exponentially more 
difficult to account for the abundance of hot clusters in the present-day universe if those same 
clusters had to be even hotter ~ 5 biUion years ago. The agreement between the measured 
abundances of hot {ksT > 6 keV) clusters and those predicted by popular cosmological 
models (sec Figures 9 and 12 in Evrard et al. 2002) would be significantly spoiled unless we 
require (AT/T) < 0.3. From equation (5) we then find r) < 0.15(rci/10 Gyr)~^7\Q^'^^. Note 
that this constraint is strongest for the hottest clusters where the total thermal energy is ~ 
]^q6'1 QYgs, and in which plausible astrophysical heating sources (such as supernovae or active 
galactic nuclei) are unable to compensate for the energy loss due to thermal conduction. Of 
particular interest is the cluster IE 0657-56 ai z — 0.296 for which the inferred emissivity- 
weighted temperature is 14.8li;2 keV, although it may be a relatively young merger (Tucker 
et al. 1998; Markevitch et al. 2002). Since heat conduction is temperature dependent, it 
distorts the shape of the abundance distribution of cluster temperatures in a generic way 
that is not degenerate with variations in cosmological parameters such as the normalization 
of the power-spectrum of density fluctuations. Inclusion of thermal conduction in future 
high-resolution hydrodynamic simulations can be used to rcflnc the above upper limit on 
f]. In these simulations, it would be important to treat the electron and ion temperatures 
separately since the thermal equilibiration timescalc is long near the virialization shock of 
clusters (Fox & Loeb 1997). 

We reiterate that the emission-weighted temperature of an X-ray cluster is dominated 
by its core which amounts to only a small fraction of the total cluster mass. A substantial 
reduction in the core temperature would merely result from the transfer of heat out of the 
hot core to the cooler, stationary cluster envelope without requiring that this heat be trans- 
ferred further to accreting material in the surrounding large-scale structure. Nevertheless, 
heat conduction is inevitably expected to preheat the surrounding gas and suppress smooth 
accretion of gas onto the gravitational potential well of the hottest clusters. If substantial, 
such a suppression would have implied that the baryon mass fraction in the hottest clusters 
is lower than its cosmic average, a result that would have been at odds with the concordance 
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model of Big Bang Nucleosynthesis and structure formation (Buries et al. 2001; Tegmark, 
Zaldarriaga, & Hamilton 2002). 

3. Conclusion 

We find that heat conduction must be suppressed by a factor rj < 0.15(rci/10 Gyr)~^Tj^^^^ 
relative to the Spitzer value, or else the cores of X-ray clusters would have cooled signifi- 
cantly over their lifetime. This suppression can be naturally produced (see Malyshkin & 
Kulsrud 2000 and §5 in Malyshkin 2001) by the magnetic fields that are inferred to exist 
in the halos of X-ray clusters (Carilli & Taylor 2002). The recent Chandra detections of 
sharp temperature jumps (cold fronts) in several clusters (Markevitch et al. 2000; Vikhlinin 
et al. 2001a) indicate an even stronger suppression of heat conduction across these jumps 
(Markevitch et al. 2000; Ettori & Fabian 2000). Similarly, a suppression factor of ~ 1-3% 
is required to account for the conditions in the interstellar gas of cluster galaxies (Vikhlinin 
et al. 2001b). If this suppression apphes also to the diffuse gas in cluster cores, then thermal 
conduction could not account for the apparent lack of cooling flows in them. Mixing of gas 
due to mergers or central heating sources, such as active galactic nuclei, could in principle 
compensate for the radiative losses in these environments (Bohringer et al. 2002; Churazov 
et al. 2002). In this paper we have shown that if the thermal conductivity had been compa- 
rable to the Spitzer value as recently suggested (Voigt et al. 2002), then cooling flows would 
have not been suppressed but rather induced in hot X-ray clusters. 

I thank Maxim Markevitch and Jerry Ostriker for useful comments on the manuscript. 
This work was supported in part by the grants 5-7768 from NASA and AST-9900877, AST- 
0071019 from NSF. 
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